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ABSTRACT 

We present a framework for the hierarchical identification and characterization of 
voids based on the Watershed Void Finder. The Hierarchical Void Finder is based 
on a generalization of the scale space of a density field invoked in order to trace the 
hierarchical nature and structure of cosmological voids. At each level of the hierarchy, 
the watershed transform is used to identify the voids at that particular scale. By 
identifying the overlapping regions between watershed basins in adjacent levels, the 
hierarchical void tree is constructed. Applications on a hierarchical Voronoi model and 
on a set of cosmological simulations illustrate its potential. 

Key words: Cosmology: large-scale Structure of Universe, methods: data analysis, 
N-body simulations, techniques: image processing 



1 INTRODUCTION 

The large scale distribution of matter observed in galaxy 
surveys and N-body computer simulations features a com- 
plex system of cell-like empty regions defi ned by a dense 



netw o rk of clusters , filaments and wal l s jKirshner et al 
19811: IColless et all 120031 : iHuchra et all 120051 ; iGott et al 



20051 ). The Cosmic Web is the result of the tidally induced 
anisotropic na ture of the gravitational collapse of density 
perturbations (|Zeldovichlll970l : iBond et al.lll996t ). 

Within this context, voids are the low density de- 
press ions from which matter is continuously draining (|lckel 
1 1984 ). Forming a key component of the Cosmic Web, voids 
emerge out of the density trough s in the primordial Gaussian 
field of density fluctuations fsee Ivan de Weygaert fc Platen! 
2009, for a recent review). As a result of their underdensity, 
voids represent a region of weaker gravity, resulting in an 
effective repulsive peculiar gravitational influence. Initially 
underdense regions expand faster than the Hubble flow and 
while they expand, matter is squeezed in between them, re- 
sulting in void boundaries consisting of sheets and filaments. 



1.1 A hierarchy of voids 

In addition to its anisotropic nature, the Cosmic Web is also 
characterized by an evident hierarchical structure. As a re- 
sult of the multiscale nature of the primordial perturbations, 
structure builds up via small scale objects into ev er larger 
structures. High resolution N-body experiments (|Springell 
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2005) display a complex and tenuous network of substruc- 
tures within the interior of voids, resembling the prominent 
Cosmic Web delineated by large haloes. The relation be- 
tween different levels in the hierarchy of the Cosmic Web 
can be defined by the voids as, at any given level of the 
hierarchy, they are the cells within which we observe the 
weblike infrastructure at the next level. 

Because of their relatively simple structure and evo- 
lution, we may better understand the gradual hierarchical 
buildup of the Cosmic Web on the basis of its void popu- 
lation. Two processes dictate the evolution of voids: their 
merging into ever larger voids as well as the collapse and 
disappearance of small ones embedded in overdense regions. 
When adjacent voids meet up and merge, the matter in be- 
tween is squeezed in thin walls and filaments, which sub- 
sequently drain towa rds the outer boundary of the voids 
l|Dubinski et all Il995f ). By identifying and assigning criti- 
cal density values to the two evolutionary void pr o cesses 
of merging and collapse, ISheth fc van de Weygaertl (|2004l ) 
managed to describe this hierarchical evolution of the void 
popul ation in terms of a two-barrier excursion set formula- 
tion (|Bond et al.lll99lh . The context of this unfolding void 
hierarchy within the Cosmic Web can be clearly understood 
withi n the Lagrangian adhesion description ( Sah ni et al.l 

[Tali. 



1.2 Reconstructing the Hierarchy of Voids 

In this study we describe our formalism for explicitly ana- 
lyzing the hierarchy of voids in the cosmic matter or galaxy 
distribution. Based on the wate rshed transform (see e.g. 
lBeucherlll982l : iPlaten et ajj|2007tl . it combines the Water- 
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shed Void Finder (WVF) with a formalism to establish the 
hierarchical structure and relationship of the detected voids. 

The g ri d-bas ed WVF method introduced by 
iPlaten et all (I2007D is able to detect voids without re- 
striction on their size and shape. A related Voronoi 
tessella tion-based imp lement at ion is the ZOBOV void 
finder (|Nevrinckl 2008). Both methods are based on the 
idea following the slope lines connecting a given point in 
space to the local minima of the valley containing that 
point. More details on the p erformance of a variet y of other 
void finders can be found in IColberg et ail (|2008h (also see 
lLavaux fc Wandeltll201oh . 

In section [5] we describe the basis of the hierarchical 
void tree formalism. The details of the technique are out- 
lined in sect. [3] We then present an illustrative test of its 
performance on a heuristic hierarchical Voronoi model in 
section [4] Its cosmological potential is outlined in sect. [5] 
followed by a short discussion in sect. H3 



2 THE HIERARCHICAL VOID TREE 

In the void hierarchy framework we identify voids indepen- 
dently at all levels of the hierarchical space, and establish 
the cross-scale relations between voids at different levels. 
For establishing the multiscale and nesting properties of the 
void ne twork, we follow the natur al path of multiscale tech- 
niques (I Aragon-Calvo et al.ll2007l ). Within this context, we 
evaluate the structure of a scalar field in N dimensions in 
an N + 1-dimensional hierarchical-space of the original field 
where the extra dimensio n represents a scale usually defined 
by a smoothing function (|liiimalll962l ; I WitkirJI 19831 ) . 

Subsequently voids between adjacent levels in the hier- 
archy are linked as a function of well-defined characteristics. 
A given parent void at the hierarchy level i is defined by 
smaller children voids at the next level, i + 1, in the hi- 
erarchy. We assign parent-child relations between voids in 
adjacent levels of the hierarchy by identifying overlapping 
volumes between the voids. A given child usually shares vol- 
ume with several parent voids higher in the hierarchy. We 
enforce a non-loop property in the hierarchical tree by as- 
signing each child void exclusively to the one parent void to 
whom the child contributes most of its volume. This con- 
straint assures that all children voids have only one single 
parent in the void tree hierarchy. 



3 RECONSTRUCTING THE VOID 
HIERARCHY 

Having established the general scheme for the void hierarchy 
tree, we need to detail its key ingredients. The first issue 
is that of the definition of the scale-space from which we 
extract the void hierarchy. The most essential element is 
the void identification at each level, which is based on the 
watershed segmentation of the scalar density - or related - 
fields. 

3.1 Scale and Hierarchical Spaces 

Proper scale spaces must have the following set of proper- 
ties: 1) linearity, 2) spatial shift invariance, 3) isotropy and 



4) causal ity. The Gaus sian filter addresses each of these con- 
straints l|Floracklll993l ). However, while the Gaussian func- 
tion is an optimal scale-space operator, it is not necessarily 
the only - or the best suited - option for the study of the 
hierarchical character of the Cosmic Web. 

The spatial filtering approach assumes that the levels 
of the hierarchy are defined purely and only on the basis of 
their corresponding spatial scale. However, it would be bet- 
ter if our definition of a a characteristic hierarchy level was 
based on the nature of the complex physical processes that 
give rise to the dark matter and galaxy distribution. Intrinsic 
hierarchical properties of the Cosmic Web such as halo mass 
functions, galaxy luminosities, galaxy morphology, etc. are 
suggestive examples. In the following, the term Hierarchical 
Spaces is used to indicate a broader class of spaces defined 
by one or more specific properties which are manifestations 
of the hierarchical nature of the Cosmic Web. This means 
they do not necessarily satisfy the requirements of a proper 
scale-space. 

In the case of N-body simulations we have access to the 
full evolution of the Cosmic Web. This allows us to control 
the relation between scales in the primordial density field. 
By using the information from the power spectrum, we can 
select those scales in the initial conditions which will grow 
and evolve faster or, alternatively, those that will not evolve 
at all. The most straightforward example would be the defi- 
nition of a linear-regime smoothing procedure that will allow 
large-scale linear fluctuations to grow while small-scale lin- 
ear fluctuations will be suppressed. This filter will act on the 
linear-regime matter distribution where all Fourier modes 
are independent and grow independently, and allow us to 
target specific hierarchy levels for further evolution towards 
collapse, ultimately producing the present-time structures. 
This low-pass filtered density field will evolve into a Universe 
with all the large scale structures in place, with their shapes 
moulded by anisotropic gravitational collapse, but lacking 
the small-scale details. 

This approach is fundamentally different from the usual 
a posteriori smoothing operation, in that it avoids the 
nonlinear effects resulting from cross talk between Fourier 
modes. It has the advantage of transparently exposing the 
hierarchy of structures imprinted in the initial density field. 



3.2 Watershed Segmentation 

The watershed transform segments an image into regions fol- 
lowing its intrinsic substructure (see lPlaten et al" I l2007l . for 
a detailed description of the method). The word watershed 
finds its origin in the analogy of the procedure with that of 
a landscape being flooded by a rising level of water: as the 
water-level rises, the watershed basins around the minima 
will ultimately meet at the ridges defined by saddle-points 
and maxima in the density field. The final result of the com- 
pletely immersed landscape is a division of the landscape 
into individual cells, separated by ridge dams. The cosmo- 
logical analogy to the landscape is suggestive: the basins rep- 
resent the underdense void regions, while their boundaries 
of sheets and ridges form the network of walls, filaments and 
clust ers that defines the Cosmic Web (|Arag6n-Calvo et al.l 
|200S|) . 
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3.2.1 Over-segmentation 

One of the practical complications of watershed segmenta- 
tion is its sensitivity to any structure, whether it is real or 
an artefact. As a result, it easily partitions a given region 
into several smaller sub-regions. This "oversegmentation" is 
commonly assumed to be the result of "noisy" structures 
superimposed on top of the more prominent - and usually 
"real" - features. In reality, the oversegmentation is not only 
set by the noise level of the image, but also by the presence 
of intrinsic and significant substructure in the field. 

The limitations of the watershed transform due to over- 
segmentation can be alleviated by the use of hierarchical 
techniques such as the hierarchical watershed l)OIsen|[l996l . 
Il997l ; lGauchHl999h . In this approach the watershed trans- 
form is computed on the image after smoothing at several 
scales or thresholding at several intensity levels. The large 
scale images will delineate large regions while smoothing 
their boundaries. The small scale images will reveal the small 
features in the image, as well as the noisy structures, while 
keeping the original boundaries. In a final step, the scale 
images are merged following a specific prescription. Often 
this involves the merging of small regions contained within 
a common parent region. 

In most watershed based hierarchical reconstruction 
schemes, the small scale images in the hierarchy are merely 
considered as an intermediate step in the reconstruction of 
the features of interest. The oversegmentation is considered 
an undesirable effect due to the noise in the image. Here we 
will use a different approach. Assu ming we may ignore the 
noise- induced oversegmentation (see lPlaten et al]|2007l ). we 
focus exclusively on the oversegmentation due to intrinsic 
structures. Instead of using only the largest scale of the hi- 
erarchy, we therefore will consider all scales simultaneously. 

3.2.2 Hierarchical Watershed 

We perform the void merging across adjacent levels in the 
hierarchy by computing only the flooding procedure on the 
watershed, i.e. without identifying the watershed bound- 
aries. This procedure segments the density field into water- 
shed basins but does not explicitly provide the boundaries 
between adjacent watershed regions. This "incomplete wa- 
tershed" focuses only on the space partitioning aspect of 
the watershed transform. This makes it straightforward to 
merge voxels between children voids on the basis of this 
incomplete watershed. This leads directly to the complete 
hierarchical void tree. After the merging procedure for com- 
pleteness we compute the full watershed transform (i.e. wa- 
tershed basins and boundaries) by performing a local flood- 
ing watershed transform restricted to the boundary voxels 
as described in lAragon- Calvo et al l l|2008h . 

Once the void hierarchy is stored in a tree structure it 
is straightforward to define functions to transverse the tree 
and extract useful information of the properties of the voids, 
their connectivity and their hierarchical relations. 



4 TEST: HIERARCHICAL VORONOI MODELS 

We tested our method with a hierarchical imple mentation of 
a Voronoi clustering model of the Cosmic Web l|Okabe!l200Cl ; 



van de Weveaert fc Icke |[l989h . This model shares similar 
spatial and hierarchical properties as the observed distribu- 
tion of matter while making it possible to objectively com- 
pare the recovered void hierarchy with the original one. Hi- 
erarchical Voronoi models used have the two main properties 
we seek to study: 1) a clear multiscale nature and 2) a hier- 
archy of nested structures. 

4.1 Implementation 

The hierarchical Voronoi model is constructed as follows: 
the top level of the void hierarchy is generated from a set 
of sparsely sampled points which define a periodic Voronoi 
tessellation. Inside each Voronoi cell we define a new set of 
points and compute the Voronoi tessellation locally on the 
points inside the cell. This local Voronoi cell is non-periodic 
and has its parent Voronoi edges as boundaries. This proce- 
dure can be repeated iteratively until the desired number of 
nested levels in the hierarchy is reached. We regularize the 
size and shape of the Voronoi cells by performing a Voronoi 
Centroid regularization on the seed points. By coupling the 
scalar and hierarchical aspects of the image we can study it 
via the canonical Gaussian scale-space. 

4.2 the Voronoi Test 

From the hierarchical Voronoi model we compute the nor- 
malized distance fie l d for each point in a regular grid (see 
lAragon-Calvo et alj l|2008l )). This field is defined as the ra- 
tio between the distance to the closest and second closest 
Voronoi seed. It yields a distance field with values of 1 at 
the cell boundaries and decreasing towards the center of the 
cell. We do this for each level in the hierarchy. Finally, all 
levels are integrated into a single distance field as X(x) as 
T(x) = y^™_ n (I' (x) /(z + 1) 2 ), where n is the number of levels 
in the hierarchy. This scale integration scheme is similar to 
the one used in othe r synthetic imag e generation algorithms 
such as Perlin noise (|Perlin. K.lll985h . By tuning the denom- 
inator of the above equation it is possible to define different 
intensity scaling relations between levels in the hierarchy. In 
our case the most prominent features in the image will be 
the largest voids. 

Next, we construct the Gaussian scale space of the im- 
age and identify voids at each scale independently. Since 
our image was constructed with two characteristic scales it 
makes no sense to use more than two smoothing scales. The 
scale-space then consists of two scales, one with no smooth- 
ing and one with a width between the size of the small and 
large Voronoi cells. 

The hierarchical merging of voids is illustrated in figure 
[T] The left top and bottom panels show the original field 
and its smoothed version respectively. The center panels 
show the corresponding watershed transform. Note that the 
smoothed field produces a distorted watershed transform. 
Both the general shape and the boundaries of the voids are 
affected by the smoothing procedure. On the other hand, 
the watershed transform of the original field reproduces the 
original boundaries between voids but it does not differenti- 
ate between levels in the hierarchy. The hierarchical merging 
of voids is shown in the top right panel. One individual void 
is highlighted in order to illustrate the void merging proce- 
dure. The parent void's area overlaps with several children 
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Figure 1. Hierarchical reconstruction of voids in a hierarchical Voronoi model. The original distance field is shown in the top left panel 
(scale 0) and its smoothed version (scale 1) in the lower left panel, the center panels show their corresponding watershed transform. An 
individual void is depicted at the largest scale in the center low panel. The hierarchical merging of the void with its children sub-voids 
is shown in the top right panel. The final hierarchical reconstruction is shown in the bottom right panel. The original shape of the large 
voids is reconstructed as well as their inner hierarchy of substructures. 



and one can see that the children that are mostly covered 
by the parent void are the ones originally inside it. The final 
result of the merging procedure is shown in the bottom right 
panel where we emphasize the large voids in the top level 
of the hierarchy (thick lines) containing smaller voids at the 
bottom of the hierarchy (thin lines). 

The hierarchical reconstruction of the voids has two 
important advantages over the single-scale watershed void 
finder: 1) It is not affected by smoothing procedures and 2) 
it explicitly gives the inner substructure of the voids. The 
reconstructed hierarchical voids contain both their original 
shape and their original level of substructure. 



We perform the linear-regime smoothing procedure 
by generating lower-resolution versions of 128 3 and 64 3 
particles from the same initial conditions. The 64 3 res- 
olution corresponds to a cut-off scale of ~ 3fa _1 Mpc, 
enough to trace voids without significant substructure. 
We followed the evolution of the box from z — 49 un- 
til the pres ent time z = 0, using the GADGET-2 N- 
body code (|Springel 2005|). From the final particle distri- 
bution we compute the density field inside a cubic grid 
of 512 voxels per dimen sion using a recent implementa- 
tion of the DTFE method (ISchaap fc van de Weygaertj2000l ; 
Ivan de Weygaert k, Schaap 1120091 '). 



5 COSMOLOGICAL APPLICATION 

We applied our algorithm to three cosmological simulations 
that are variants of the cold dark matter scenario. The sim- 
ulations cover the three possible geometries of the Universe: 
flat, open and closed, with cosmological parameter of (fi m , 
n A ) = (0.3,0.7) for the the flat ACDM model, (0.1,0.7) for 
the open ACDM and (0.5,07) for the closed ACDM Universe. 
Each simulation consists of 256 3 dark matter particles in a 
200fa -1 Mpc box. All simulations sh are the same Hubb le pa- 
rameter, h = 0.7 and as — 0.8 (see lArava-M clo 2008, for a 
detailed description). 



The size distribution of voids in different cosmologies 
and levels of the hierarchy are shown in figure [2] The mean 
void sizes of voids in all cosmologies are 11,13 and 18 Mpc/ft 
for levels 0,1 and 2 respectively. While the voids at the top of 
the hierarchy (level 2) are clearly the largest, the mean size 
and distribution of voids in levels and 1 are very similar. 
All distributions in the three cosmologies have similar peaks 
at a given level in the hierarchy. However, there are differ- 
ences in the overall shape of the distributions. Compared to 
the LCDM and SCDM, the OCDM universe has a higher 
tail towards large voids. The fact that the order OCDM- 
LCDM-SCDM is observed in all the distributions gives us a 
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Level Level 1 Level 2 




Figure 2. Distribution of void sizes for three cosmologies: SCDM (solid) ACDM (dashed) and OCDM (dotted) computed at three 
different levels of the hierarchy going from the (left) bottom of the hierarchy (smallest scale) to (right) the top of the hierarchy. 



good indication of the ability of our method to discriminate 
between cosmologies. 



6 CONCLUSION AND FUTURE WORK 

We introduced a framework for the identification of voids 
and their hierarchical properties. The hierarchical nature 
of the void network makes our method a powerful tool for 
its description and characterization. The Hierarchical Void 
Finder shares the advantages of the Watershed Void Finder, 
while addressing some of its limitations such as the over- 
segmentation and the reconstruction of the void boundaries 
after strong smoothing of the density field. 

In order to test our method we introduced a hierarchical 
implementation of Voronoi models. These heuristic models 
share the multiscale, hierarchical and topological properties 
of the Cosmic Web. As such the hierarchical Voronoi mod- 
els represent a valuable tool for testing algorithms for LSS 
analysis. 

We extend the idea of scale-space in order to account for 
non linearities and physical processes. We discuss a Gaussian 
smoothing in the initial conditions. By applying, a Gaus- 
sian smoothing in the linear regime before there is cross- 
talk between Fourier modes we are able to cleanly expose 
the hierarchy of structures in the evolved non-linear matter 
distribution. 

In a following paper we will describe the basis of the 
hierarchical space and explore in more detail the properties 
of the void network and its potential for constraining cos- 
mological parameters. 
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